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, We calculate numerically scattering phases for elastic meson-meson 

scattering processes in the strongly coupled massive Schwinger-model 
with an SU(2) flavour symmetry. These calculations are based on Lu- 
scher's method in which finite size effects in two-particle energies are 
[ exploited. The results from Monte-Carlo simulations with staggered 

■ fermions for the lightest meson ("pion") are in good agreement with 

the analytical strong-coupling prediction. Furthermore, the mass spec- 
trum of low-lying mesonic states is investigated numerically. We find 
a surprisingly rich spectrum in the mass region [m T ,4mJ. 

Q-i, 

Jh ; PACS: ll.10.Kk, 11.15.Ha 

Keywords: Lattice gauge theory; Schwinger model; mass spectrum; 
elastic scattering phases 



X 

1 Introduction 



Present experiments at high energy accelerators are mainly scattering ex- 
periments. A task for theoretical QCD is to calculate quantities which can 
be compared to the data of such scattering experiments. In deep inelastic 
scattering processes for example, reliable results are obtained by perturba- 
tive calculations. But problems occur in experiments with small momentum 
transfers. These are, e.g., resonant elastic scattering processes, like the 
occurrence of the p-resonance and the A-resonance in it-it and n-p respec- 
tively scattering processes. For such low-energy phenomena perturbative 
calculations in the strongly-coupled QCD region are not suitable. 
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By contrast there exists an appropriate method in the framework of lat- 
tice calculations proposed by Liischer in which scattering phases of the con- 
tinuum theory can be calculated for elastic scattering processes in massive 
quantum field theories . In this method one makes use of the fact that 
for large but finite spatial extension the volume dependence of the energies 
of two-particle states is determined by the S-matrix for elastic scattering 
processes in infinite volume. As the determination of energies in finite vol- 
umes is possible in Monte-Carlo simulations this method is very useful for 
lattice calculations. 

There are many successful applications of this method for the calcu- 
lation of scattering phases in bosonic models in two and four dimensions 
[PHlO|]. The basic formulae for Liischer's method have also been derived 
in a fermionic theory fLl| . The agreement of the numerical data with the 
predicted scattering phases for the fermion-fermion scattering in the Gross- 
Neveu model in two dimensions [12| confirms the usefulness of Liischer's 
method also in fermionic models. 

Our aim is the application of Liischer's method to the determination of 
elastic scattering phases in a meson-meson system of the massive Schwinger 
model with an SU(2)j flavour symmetry in the continuum. The experience 
gained in this project shall support future investigations in the QCD, e.g. of 
the p-resonance in the ir + ir~ — > p° — > ir + ir~ scattering process. In this con- 
text the Schwinger model has some useful features: The Schwinger model is 
a fermionic gauge field theory which has many properties in common with 
QCD, such as confinement, the U(1)a anomaly and a non-trivial vacuum 
structure. In particular the massive Schwinger model with Nf = 2 offers a 
complex mass spectrum and therefore conditions (and problems) comparable 
to QCD with ii,d-quarks as far as the mesonic energy spectrum is concerned. 
Furthermore, the simulations for the determination of the scattering phases 
are extensive, e.g. because of the calculation of fermionic eight-point func- 
tions. Therefore a two-dimensional model has the advantage to be not as 
expensive to simulate as the four-dimensional QCD. 

Another advantage of the massive Schwinger model with Nf = 2 is that 
there exist analytical predictions for strong couplings: The massive SU(2)j 
Schwinger model is not analytically solved. But for strong coupling the "pion 
sector" of this model, i.e. the sector in which the lightest meson-triplet 
("pions") occurs, can be approximated by the sine-Gordon model. In the 
sine-Gordon model the mass spectrum and the elastic S-matrix have been 
calculated ||l3|-[l5|| . The primary aim of our investigations is the comparison 
of the predictions for the scattering phases from the sine-Gordon model with 
the numerical results in the massive SU(2) ^ Schwinger model calculated with 
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the procedure of Liischer. By this comparison we also check ambiguities, the 
CDD-poles [|l5|], in the analytical result for the S-matrix. 

For a successful determination of the scattering phases it is necessary 
to have a good knowledge of the mass spectrum in the massive SU(2)j 
Schwinger model: 

First, investigating the particle masses in the pion sector we determine 
the strong coupling region in which the approximation of the pion sector by 
the sine-Gordon model and thus the prediction for the scattering phases of 
the pion-pion scattering are nearly exact. Secondly, the analysis of the mass 
spectrum is useful to estimate the energies of possible additional scattering 
processes besides the it-it scattering and of inelastic thresholds in the energy 
region of the elastic it— it scattering. Thirdly, we investigate systematically 
lattice artifacts — finite size and 0{a) effects — affecting the particle masses. 
Such investigations are very important because small systematical deviations 
in the two-particle energies lead to large errors in the scattering phases 
calculated from these energies. 

The paper is organized as follows: In section || the known analytical 
results for the mass spectrum and the scattering phases are summarized. 
In section ^ the numerical methods, especially an improved method for the 
determination of energies from correlation matrices, are discussed. The rela- 
tion of the symmetry groups on the lattice and in the continuum is described 
in section |||. In section ||| the numerical results for the mass spectrum of the 
massive SU(2)j Schwinger model are presented. The numerical results for 
the scattering phases for the elastic it-it scattering are compared to the 
analytical prediction in section ||. In this section also the contributions of 
single-particle states and of states of two non-interacting pions to four- 
meson correlation functions are discussed. The paper ends with a summary 
of the main results. 

2 Analytical predictions in the massive SU(2)^ Schwinger 
model 

2.1 Bosonization 

In Minkowskian space-time the action of the massive Schwinger model with 
a SU(2)^ flavour symmetry is 1 : 

S = f d 2 x\^-- A F, v F^ + ^(t?-e4-mo)^j, (1) 

1 H = c = 1, g 00 = 1, gn = -1. 
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This model has many properties in common with the four-dimensional QCD: 
In the massless case (mo = 0) there exists an anomalous axial U (l)-current 
[ |T(| and a massive meson ("77" -meson) which corresponds to the 7/ -meson 
in QCD. Also the vacuum structure of the Schwinger model is not trivial. 
Analogously to the screening of colour charges of test particles in QCD 
with dynamical quarks the electric charge of test particles in the Schwinger 
model is completely screened. Furthermore, the "quarks" which denote the 
fundamental fermions in (|l]) are confined in the massive Schwinger model 
with Nf = 2 for |0| 7^ ir/2. Thus there exist no free states with a non-zero 
fermion number |17fl . 

The basis for the bosonization of the Schwinger model 2 is the equivalence 
between the sine-Gordon model and the massive Thirring-model |l8|]. This 
equivalence gives the possibility to express the fermionic fields ip 1 , if) 2 (and 
also the gauge fields by means of the equation of motion) in (|l]) by two 
bosonic fields 0+ and <f>_. The bosonized form of (|l]) which is obtained this 
way was derived in ref. fl9| . In Hamiltonian formulation one obtains: 



2 > 2 + ^U+^+^ 12 
-2cm^/ 2 fi 1 ^ 2 : cos (VZk^+j : m : cos ( v / 27r<^- 



■mo 



: inmo denotes normal ordering with respect to a free field with mass \x and 
mo respectively. n + and II_ are the conjugate momenta of the Bose fields 
4> + and <fi_. The parameter /i is defined as fi = e^j2jix. The constant c 
stems from the bosonization formulae of the massive Thirring model and 
has the value c = exp(7)/(27r) (2(| (7 is the Euler constant). 

We now confine ourselves to the strong coupling limit of the model, i.e. 
to the case where mo/e <C 1. Expanding the cosine-terms in (|2|) the mass 
scale of the (/>_- and </> + -fields is given by (m^ 2 /j. 1 ^ 2 ) 1 ^ 2 and \i respectively. 
For mo/fi <C 1 the (/>+-field has a much higher mass in (0) than the 0_-field. 
Hence for mo/e small the </>+-field has only small effects on the dynamics 
of the </>_-fields so that the interaction terms in (|2|) are negligible 

if only the </>_-field is considered. With these approximations one finally 



obtains with small mo/e the following theory for the </>_-field [19]: 

^_ = : ^n 2 _ + i(d 1( M 2 -^cos(/?5G0-) W, (3) 



2 In the following the term Schwinger model denotes the massive case with an SU(2)^ 
flavour symmetry. 



4 



m' = 2^c 2 / 3 (^) ' e, (3sg = V^. 

This Hamiltonian represents the sine-Gordon model with a coupling param- 
eter Psg- Hence it is possible to make use of the analytical results for the 
sine-Gordon model in order to obtain predictions for the lightest particles 
in the Schwinger model. The approximation of the Schwinger model by the 
sine-Gordon model is, of course, only valid in the c/>_-sector ("pion sector") 
of the Schwinger model for mo/e small. 



2.2 Mass spectrum and scattering phases 

The sine-Gordon model is an integrable model. Solutions for the mass 
spectrum in the quantized theory have been first derived in a semi-classical 
approximation The mass of the soliton (A) and antisoliton (A) is 

predicted to be 

8m ' 111 
M A = M A = - ? =-^- S . (4) 

The number of additional stable particles B n , n = 1,2... < 871-/7' depends 
on the parameter fisG- In the case Psg = there exist two stable states 
B\ and B2 with masses M\ = Ma and M2 = Hence for mo/e — > 

there are four particles in the Schwinger-model which are related to the </j_- 
sector. The SU(2)j flavour symmetry of the Schwinger model is visible in 
the sine-Gordon model for (3sg = V / 2vr as the particle spectrum consists of a 
triplet {A,A,B\) and a singlet (i?2)- The quantum numbers (Isospin, Parity, 
G-parity) of the triplet are 

jPG _ 1 -+ and of the singlet jPG _ ++ gjj 

In analogy to QCD we call the pseudoscalar triplet "pion" and the ++ - 
particle "/o" -meson. From (g) and @ one obtains for the pion mass in the 
Schwinger model 

= i066(^) 2/ % (5) 



and for the /o-mass 



. (6) 



Exact calculations for the mass gap in the sine-Gordon model in the full 
quantum field theory show that the semiclassical results are nearly exact 
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[ p~f . p| ] . The difference of the exact result 



rrv ~ 2.008 ^) 2/3 e (7) 

for the pion mass to the semiclassical formula (g) is about 3%. 

In the original derivation of the Hamiltonian (p|) in ref. |l9| the param- 
eter m' contains the additional factor (cos(0/2)) 2 ' 3 with G G [— it, ?r]. This 
factor reduces the pion mass in (|5|) and ([?]) for 0^0. All our Monte-Carlo 
simulations yield numerical data for the pion mass which are above the val- 
ues of (||), indicating ~ 0. Therefore we put = in all relevant formulae 
in this paper. 

Besides the particles which have their origin in the 0_ -sector there exists 
a particle in the Schwinger model resulting from the ^ + -field in (|2|). It 
has quantum numbers I PG = ("r/"-meson) pll. In the massless limit 



mo — ► it becomes the massive boson of the massless model with mass [17] 



/2 

m = = ey — • (8) 

In the massive theory we therefore expect for mo/e <C 1 a mass of the 
order m v ~ \i. Other states whose mass is above the mass of the n-, fo~ 
and 7y-meson can exist because of the interaction and the <^>+ self- 

interactions in Many of such possibly unstable states have been found 
in our numerical simulations (see section |5.4| ). 

To obtain an analytical formula for the elastic scattering phases for the 
pion-pion scattering in the strongly coupled Schwinger model one can once 
again make use of its approximation by the sine-Gordon model. The S— 
matrix in the sine-Gordon model has been derived exactly [15|. The calcula- 



tion is based on the unitarity and crossing-symmetry as well as on the factor- 
ization property of the S-matrix in the sine-Gordon model. For (3sG = V^n 
one obtains the same scattering matrix element S(9) for the three possi- 
ble scattering processes with soliton (A) and antisoliton (A): AA — ► AA, 
AA — * AA, AA — > AA: 

sinh(fl) + i sin tt/3 

b[ V' = ■ \Ta\ ~ Ti ■ ^ 

smn(t7j — i sin %/ 6 

The difference of the rapidities 6\ and 02 of the two particles is denoted by 9: 
9 = 9\ — #2- It is related to the particle momentum k by 9 = 2 arsinh^/M^). 
Equation @ is an approximation for the scattering matrix element for the 
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elastic pion-pion scattering in the strongly coupled Schwinger model. Defin- 
ing scattering phases 5(6) via 



S(Q) = e 2i5 ^ (10) 

one obtains the connection 5(k/m n ) between the elastic scattering phases 
and the particle momenta. The aim of our numerical investigations is to 
compare the scattering phases 5(k/m,T T ) obtained by the method of Liischer 
with Monte-Carlo simulations with the analytical results @. 

The scattering matrix element @ is unique up to multiplicative "CDD"- 
terms: 

£± smh{6) — i sm(afc) 

Expression (^) is the "minimal" solution with the lowest number of poles. 
Arguments that this minimal solution for the soliton-soliton scattering in 
the sine-Gordon model is exact were given in ref. There exist two 

states in the sine-Gordon model for (3sg = v / 27r with masses Ma and \[?>Ma 
which can be interpreted as bound states of soliton and antisoliton. This 
corresponds to the two poles 6 = iir/3 and 6 = i2n/3 in the scattering 
matrix element (|9|) of the soliton-antisoliton scattering. Furthermore, the 
minimal solution is in agreement with the results for the S-matrix in the 
semiclassical limit fisc ~ > 0- Our numerical results discussed in section [Jl] 
also support the scattering phases of the minimal solution. 

Because of the factorization property of the S-matrix the scattering ma- 
trix elements of scattering processes involving the states B n can be calcu- 
lated by using the scattering matrix element of the soliton-soliton scattering. 
The states B n can be interpreted as bound states of soliton and antisoliton. 
Hence arbitrary scattering processes can be understood as a sequence of scat- 
tering processes between soliton and antisoliton. Calculating for fisc = V%ft 
the scattering matrix elements for the processes B\A — > B\A, B\A — > B\A 
and B\B\ — > B\B\ one obtains in each case the expression (^). Hence all 
scattering processes within the triplet (A, A, B\) have the same scattering 
matrix element. That means that the elastic scattering phases for the tt-tt- 
scattering in the Schwinger model for m^fe — > are identical for all isospin 
channels. 
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3 Numerical methods 



Our starting point for Monte-Carlo simulations of the Schwinger model is 
its Euclidean lattice action with Kogut-Susskind fermions: 

S = S w + Y,XxM{x,y) X y, (12) 

x,y 

Sw = /3^{1-Ref/ P (x)} , 

X 

U p (x) :=U^{x)Ut{x + 2)U 2 {x + l)Ui(x) , 
M (x,y) = -^^2vfi(x){U{x,x + fi)5 x+fl)y -U{x,x - fi)S x ^ y } 

+m a 2 6 X:y , mix) = 1, 7] 2 {x) = {-l) Xl . (13) 

The sum over the space-time coordinates x = (xi,x 2 ) runs for the space 
coordinate X\ = x from a to La and for the time coordinate x 2 = t from a 
to Ta. We use the Wilson action with the link variables U^(x) as a compact 
formulation of the gauge fields: U^(x) = U(x + /t, x) G U(l). The Kogut- 
Susskind fermions x, % on the lattice have one flavour so that a two-flavour 
Schwinger model in the continuum is simulated. In the simulations we use 
periodic boundary conditions for the link variables and periodic boundary 
conditions in space and antiperiodic in time for the fermionic fields. 

In the naive continuum limit the parameter (3 has to be chosen (3 = 
l/(e 2 o 2 ) to obtain the continuum action of the Schwinger model. Therefore 
for finite a we define a coupling e on the lattice by e := /? -1 / 2 . Here and in the 
following the lattice constant a is set to one. Dimensional and dimensionless 
parameters and fields are denoted by the same symbols provided there is no 
confusion. 

3.1 HMC and topological ergodicity 

We use the Hybrid Monte Carlo (HMC) method with pseudofermions and 
the conjugate gradient algorithm for the inversion of the fermion matrix. 
The HMC provides an efficient algorithm for the calculation of Euclidean 
path integrals in the Schwinger model. For (3 < 5 it appears to be ergodic 
with respect to the transition between different topological sectors ("topo- 
logical ergodicity"). Topological sectors are characterized by the topological 
charge Q of the gauge field configurations. We use the geometrical definition 
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P=5.0 

m =0.15 

32x32 



p,= 10.9(1.1)% 
N=7204 



Figure 1: Probability distribution and time series of the topological charge Q for 
a typical simulation with N — 7204 configurations for (/3 = 5, mo = 0.15). 
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Op S (— 7r,7r] is the plaquette angle which is defined by U p (x) = exp(i# p ). 
The tunneling probability pr, he. the probability for a change of the topo- 
logical sector in the HMC simulations, is about 10% for (3 = 5 and for ten 
trajectories between the measurements, reaching an acceptance probabil- 
ity of about 60-70%. The probability distribution of the topological charge 
p(Q) for this coupling is plotted in fig. [I]. Its shape can be understood 
phenomenologically from investigations in the pure U(l) theory p4]. The 
mean action (and also the minimum action) in each topological sector is 
proportional to Q 2 on large lattices. Assuming a Boltzmann distribution 
one obtains a Gaussian like shape for p(Q) as in fig. |l]. The difference of 
the mean (and minimum) value of the action for adjacent topological sectors 
increases for increasing (3. Therefore the width of p(Q) becomes smaller for 
larger (3 values so that for /3 > 10 almost only the sector Q = contributes. 

In the quasi micro canonical calculation of the trajectories in the Hybrid 
Monte Carlo the action serves as a potential. Hence one is confronted with 
the problem that for high /3-values the tunneling between adjacent topolog- 
ical sectors is suppressed. In fact the tunneling probability px for (3 > 10 is 
far below 0.1%. In that case the HMC algorithm is not suitable to perform 
simulations having topological ergodicity, which is also known from QCD. 
This is, however, only a practical problem for the calculation of ensembles 
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of moderate size which becomes less relevant on very large ensembles. 

Another problem follows from the fact that the correlation functions have 
different values in different topological sectors: In contrast to investigations 
for (3 = 1 |^] it turns out that for high values of (3 the pion mass differs 
significantly for different topological sectors p6|| . Also the exponential decay 
of the correlation function C(t) for the pion is modified for Q / (see 
ref. (2?]] for Wilson-fermions) . For large t and \Q\ the values of C(t) are 
substantially smaller than it is expected for a pure exponential decay. Thus 
it is only possible for small values of t and \Q\ to determine the pion mass by 
fits with exponential functions. But one is interested in high /3-values as the 
continuum limit a — > of the Schwinger model requires (3 — > oo. A possible 
solution to the problems is the introduction of improved HMC algorithms 



which force the tunneling between different topological sectors |28, 29[. 

We have chosen an alternative way by calculating expectation values 
for (3 > 10 only for Q = 0. We started the Monte-Carlo simulations with a 
suitable gauge configuration with Q = 0. Because of the small px for j3 > 10 
the generation of ensembles of 0(50000) was possible without tunneling 
into sectors Q / 0. This procedure is justified because the probability 
distribution p(Q) of the topological charge for (3 > 10 is quite narrow. Hence 
the expectation value of an operator O which can be decomposed on the 
lattice as 

oo 

(o)= (o %/],*> ( 15 ) 

u=— oo 

is dominated by the Q = sector: 

(O) ~ (O5 Q[ul0 ) . (16) 

This approximation is valid as long as the absolute size of the expectation 
values for Q ^ in ([15]) is much less than for Q = 0. A possible reason for 
large expectation values for Q ^ could be the influence of approximate 
zero modes of the fermion matrix M on fermionic correlation functions like 

(XxXxXyXy) = (M- 1 (x,x)M- 1 (y,y) - M- 1 (x,y)M- 1 (y,x)) u , 

{...) u = ^ J VUdetM ...e~ Sw , 



VUdetMe 



-Sw 



with the consequence that the elements of the inverse fermion matrix become 
large. But as the fermion matrix for uiq = is antihermitian for staggered 
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fermions it has no zero modes for rriQ ^ 0. Therefore we do not expect 
large contributions for Q 7^ for the bare masses rriQ / we used in our 
simulations. (The same problem for the massless Schwinger model with 
staggered fermions was discussed in ref. With the approximation (|l~6| ) 

the quantity ( O ) can be calculated in two ways: Either the expectation 
value of O 8q\u\q for all possible gauge field configurations is determined or 
the calculation of O is performed with the probability density 

P Q= o = ^S Q[u]t0 e- s . (17) 

In the latter case the acceptance probability for configurations with Q[U] ^ 
is zero. Thus only the ergodicity of the algorithm within the topological 
sector Q = is needed. In real simulations the explicit rejection of a con- 
figuration according to the <5-function in ([l7]) is not necessary for j3 > 10. 
As it was mentioned above the tunneling probability pr for high (5 is so low 
that no tunneling occurs for ensembles of moderate size. 

A known disadvantage of HMC simulations for high /3-values is that due 
to weakly fluctuating gauge fields the autocorrelation of gauge and fermionic 
operators increases. Therefore for (3 = 10 we have calculated ten HMC 
trajectories between the measurements to obtain uncorrelated values for the 
interesting operators. 

3.2 Enlargement of correlation matrix 

The main task of the analysis of the numerical data is the extraction of 
energies from connected correlation matrices Cij(t) = (Oi(t)Oj(0)) c , i,j = 
1 . . . r. A standard method for this is the calculation of generalized eigen- 
values of the matrix C(t) = (Cy(t)): 

C(t)w^(t) = AKt)C(t )^ (0 (i) , I = 1 • • -r , (18) 

to <t , 

The generalized eigenvalues A/(t) have the following form for large times 
(t-oo)[@: 

Aj(f) = (tr,)*-*^-^-*)) + ( e -AB,(*-*o)) , Gl = ±1 . (19) 

AEi is the smallest difference between the energy Ei and other energies in 
the spectrum: AEi = min^ \Ej — E\\ . The prefactor (d;)* - * is due to 
the usage of staggered fermions. If o\ = +1, —1 we call the eigenvalues and 
energies non-alternating and alternating respectively. 
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Figure 2: Logarithmic plot of the generalized eigenvalues of a four-meson corre- 
lation matrix Cij(t), i,j = 1 . . .4. Positive values are represented by circles and 
negative values by crosses. 



One possibility to relate the numerically found eigenvalues \(t) for differ- 
ent times t to a fixed I is to consider the absolute size of the eigenvalues. This 
may lead to problems if the energies in the spectrum and hence the eigenval- 
ues are close to each other. A procedure which turns out to be often more 
successful is the assignment by means of the generalized eigenvectors w^: 
The generalized eigenvectors w^ l \t') of a timeslice t 1 (normally t' = to + 1) 
are reference vectors for the assignment of the eigenvalues. A generalized 
eigenvalue Awj^t) on the time slice t is assigned to the eigenvalue Aj(t') if 
7T G S r is a permutation which fulfills: 



E 

i=i 



IV 



(0(t').^«0)( t ) maximal . (20) 



This means that the eigenvectors of the time slice t should be chosen as 
parallel as possible to the reference eigenvectors of the time slice t'. 

In some cases even the above methods are not sufficient to extract ener- 
gies from a correlation matrix. In figure || the generalized eigenvalues for a 
typical 4x4 four-meson correlation matrix are plotted. It is apparent that 
the different eigenvalues are close to each other and that they have time os- 
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dilating corrections which are due to alternating energies. This could lead to 
a wrong assignment of the eigenvalues. Apart from the highest eigenvalue 
it is also not possible to perform a fit with a single exponential function. 
Analysis also for other parameter values shows that there are only slight 
improvements if the number of operators 0j(t) is increased. 

In such cases a substantial progress can be made by enlarging the corre- 
lation matrices. In this (to our knowledge) new method a correlation matrix 
is enlarged to its double size to reduce the contributions of higher energies: 

Let C(t) be a (not necessarily symmetric) r x r correlation matrix with 
elements Cij(t) of the form: 



2r 
k=l 



(21) 



A (2r) x (2r)-matrix E(t) is now constructed by 



E. 



2i-l+i,2i-l+j 



Cij(t + i + j) 
i,j = l...r, 
U = 0,1 . 



(22) 



E.g., a 4 x 4 matrix E{t) which is made up of the elements of a 2 x 2 matrix 
C(t) has the following form: 





( Cii 




Cn(t + 


1) 


Ci 2 (t) 


Ci 2 (t 


+ 


1} \ 


E(t) = 


Cii 


(t + i; 




2) 


C 12 (t+1) 


C 12 (t 


+ 


2) 




(*) 


C 2 i(t + 


1) 


c 22 (t) 


C 2 2(t 


+ 


1) 




V c 21 


(t + i; 


C 2 i(i + 


2) 


c 22 (t + i) 


C 2 2(t 


+ 


2) / 


Defining vectors 


( e '(fc) 


,e«) 


via 












E i iji(t)= 


=£« 


'(fe)Jfc) 






,i' = l 




2r , 



(23) 



k=l 



one obtains for these vectors with definition d2 



„'(*) 

-2i-l 

„(*) 
-2i-l 



/(fc) 

V > 

„(fc) 



"2i 
(fe) 

4- 



(24) 



The form ( p3[ ) of the elements of the matrix E(t) is suitable to determine the 
energies Ej~ , k = 1 . . . 2r by the calculation of the generalized eigenvectors 
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Figure 3: Generalized eigenvalues of an enlarged 4x4 correlation matrix. Plotted 
are those eigenvalues which correspond to the ones in fig. ^ as well as an additional 
alternating eigenvalue. 

of E(t). For these calculations it is required that the vectors e ,<yk \ e^ k ' are 
linearly independent. This has to be checked from case to case during the 
numerical analysis. 

The main advantage of the enlargement of the correlation matrices is 
that according to ( |24D the amplitudes e - , , j = 2, 4, . . . , 2r are weighted 
with the factor exp(— This results in a suppression of higher energies in 
the spectrum by a stronger exponential factor in (|23|). Hence one obtains a 
smaller distortion of lower eigenvalues through higher energies and a better 
separation of the generalized eigenvalues. Furthermore, on all time slices 
a better differentiation between alternating {at = —1) and non-alternating 
(<7fc = +1) eigenvalues is possible as the eigenvectors change qualitatively 
for <7fc = — 1 because of the multiplication of each second element with a^. 
The success of this method is obvious in fig. || which is the analogue of fig. 
^ for an enlarged correlation matrix. An improvement of the determination 
of the eigenvalues and energies is obvious. A further effect is the doubling 
of the number of determinable energies for a correlation matrix with fixed 
size. 

In real applications the time extension T of the lattices is finite and 
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a replacement of exponential functions by hyperbolic function in (^T]) is 
necessary. The adjustment of the enlargement method in this case is done 
by a modification of the equations (|22|): 



E2i-l,2j-l(t) 
E2i,2j-l(t) 
E2i-l,2j{t) 
E 2i,2j(t) 




(25) 



and (H): 




/(*) 




a k cosh(-E k )v- 
a k cosh(-E k )v\ 



(fc) 



(fc) 



(26) 



Five matrices C(i — 2), C(t — 1), . . . , C(t + 2) are needed per constructed 
matrix E(t). Therefore the number of time slices on which the matrix E(t) 
can be evaluated is reduced by four. Concerning the matrix enlargement one 
also has to be aware of the fact that in comparison to the matrix C(t) larger 
statistical errors for the generalized eigenvalues of E(t) have to be taken into 
account. The reason is that the signal/noise-ratio of a general correlation 
matrix C(t) decreases for larger times t. As the correlation matrix C(t + 2) 
with its statistical error contributes to the matrix E(t) the signal to noise 
ratio of E(t) is in many cases worse than that of the matrix C(t). This is 
true as far as the matrices C(t) are statistically uncorrelated for different 
time slices t. If they are statistically correlated the opposite effect might 
occur. In that case the calculation of the generalized eigenvalues can result 
in a cancellation of the statistical fluctuations of the matrices C(t). This 
might lead to an underestimation of the statistical error. Therefore a careful 
check of the statistical errors is necessary when applying the enlargement 
method. 

4 Symmetry transformations 

In all simulations we used time slice operators in which the dependence on 
the space coordinate drops out by summation. Therefore we are interested 
in those symmetry transformations of the action which leave the time t 
invariant. We use in particular the "bosonic" representation D, which is 
defined by its action on gauge-invariant, mesonic operators — composed of 
fermion, antifermion and gauge field — with total momentum P = 0. In this 
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representation transformations that multiply fermion and antifermion with 
phase factors exp(— i<f>) and exp(z^) respectively are trivial. Such symmetries 
which are not trivial in this representation form the continuum time slice 
group CTS and the lattice time slice group LTS respectively. In the following 
a short description of the groups in two dimensions and the connection 
between lattice and continuum irreducible representations is given. The 
procedure to connect the lattice and continuum groups is presented in more 
details in refs. [0-0. 

The relevant symmetries of the continuum action of the Schwinger model 

are 

• parity 3?, 

• G-parity^ := e^ 2 ^ 

= io-$> , "rf charge conjugation, 



• SU(2)j flavour transformation. 

The G-parity is defined as charge conjugation with subsequent rotation 
in isospin space. The above symmetry transformations for the contin- 
uum action commute. They are the generators of the free group CTS = 



^^SU(2)j ). Following the notation of ref. [35] we characterize an irre- 



ducible bosonic representation (also called symmetry sector) of the CTS by 
A^ pCTG . The quantum numbers of the discrete transformations 3 d and £f are 
denoted by ap,aa = ±1 and the irreducible representation of the SU(2)j 
group by D. 

For the Schwinger model on the lattice with staggered fermions the fol- 
lowing symmetry transformations exist: 

• shift in space-direction 

• spatial inversion <y, 

• charge conjugation ^o- 

The exact action of these symmetry transformations on the fermions and 



gauge fields is described in ref. |3l|]. According to the continuum case we 
denote an irreducible bosonic representation of the group LTS = ( J* ) 
by 

&n*i<rc = ( ai ,a I: ac) ■ (27) 
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In the following mainly the short notation on the right-hand side of equation 
(§7|) is used. 

The ansatz for the embedding of the symmetry group of the LTS in the 
CTS are the equations 

^(2x) = 2^/ 2 Y,lf f ) U{2x,2x + H)x 2x+ Hi (28) 

H 

^(2x) = 2-V 2 Y,X 2x+ HU(2x + H,2x) 1 f f Y , (29) 

H 

# = 0,1,2,12 , £ = 0,1,2,1 + 2 , 1H = 1,71,72,7172 , 
which connect the fermionic fields on the lattice x, x with those of the con- 



tinuum if),ip [36]. By applying the symmetry transformations of the LTS 
on the fields if) in terms of the lattice field \ one obtains transformation 
properties of the continuum fields. From this one obtains the connection 
between the CTS and the LTS. 

For the determination of the mass spectrum we use mesonic operators, 
i.e. tensors of rank two 

and for the scattering phases tensors of rank four 



ip f fl = ip f ip f ' (30) 



4] = ^ty W . (31) 

It is known that the irreducible representations of the CTS restricted on the 
subgroup LTS are reducible and decompose into a direct sum of irreducible 
representations of the LTS. For the tensors (30) and ( |3l| ) the connection 
between the irreducible representations is listed in table [I] which is the cor- 
rected version of table 1 in ref. [25]. Each lattice sector is connected with 
two and four respectively different continuum sectors. 

A better assignment of an energy determined in a definite lattice sector 
to a sector in the continuum is obtained by using the results of the transfer 
matrix formalism for staggered fermions (see ref. [37| and the references 



therein). With these results the factors ±(±1)' in the eigenvalues of the 
correlation matrices C(t) can be related to discrete quantum numbers of 
the LTS and CTS. Hence for tensors of rank two a unique assignment of 
the energies obtained for a definite lattice sector to a state in a continuum 
sector is possible. 

In contrast to the singlet states there exist three lattice sectors for tensors 
of rank two in which a triplet with given continuum quantum numbers can 
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LTS 


CTS 


fl <7/ 


D ap oq 


'rank 9^ 

D ap oq 


D ap a G 


A ff P ff G 

S 

-D <7p 0"G 


yon Lr /1l 

D ap ac 


D ap uq 


+ 

- + - 
+ + 

+ 

+ + 
+ + 
+ + + 


1 + 
1 + + 
1 

1 + 
1 + 
1 + + 
1 

1 + 


1 

1 + 
1 + 
1 + + 


+ 
+ 
+ + 


2 - - 
2 + 
2 

2 + 
2 

2 + 
2 

2 + 


2 + 
2 + + 
2 + 
2 + + 
2 + 
2 + + 
2 + 
2 + + 


1 

1 + 
1 

1 + 
1 + 
1 + + 
1 

1 + 


1 + 
1 + + 
1 + 
1 + + 


+ 
+ 
+ + 



Table 1: Connection between the irreducible representations of CTS and LTS. 



be investigated. Therefore it is much more easier to find a suitable operator 
for the investigation of the triplet states, e.g. the pion. According to 

1"+ ® 1 _ + = 2++ © 1 ++ © ++ 

a two-pion state exists in three different isospin channels. This leads to four 
different lattice sectors in which a priori investigations for the it-it scattering 
can be performed. As it will be shown in the discussion of the numerically 
determined scattering phases in section ^ there exists a posteriori only one 
lattice sector which is suitable for the determination of the elastic tt-tt scat- 
tering phases. 



5 Mass spectrum 

In the context of this paper the mass spectrum of the Schwinger model is 
important for two reasons. First, the lightest particles offer the opportu- 
nity to check the analytical predictions (in the strong coupling limit) for the 
mass spectrum in the Schwinger model with numerical methods. Secondly, 
no exact calculation for the complete mass spectrum of the massive Schwin- 
ger model exists up to now. But for the interpretation of the numerically 
determined scattering phases it is necessary to have information about the 
lightest particles in the spectrum. Particles with masses m > 2m n might 
occur as resonances in the elastic tt-tt scattering whereas other mesons with 
m < 2???,^ might cause inelastic processes at energies below the threshold 
E = Am-n. 
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5.1 Correlation functions 



The simplest time slice operator which can be used for the calculation of 
the pion mass is 

= jE^-ir, (-ir = (-i) x+t • (32) 

X 

The easiest way to calculate the correlation function C(t) = (O(t)O(0)} for 
the investigation of a triplet is to consider only the connected part in the 
correlation function: 

C^^i^dM-^x,^)! 2 ), . (33) 

It is expected that for exact flavour symmetry in the continuum limit, the 
CV-function yields the correct pion energies [38|. 

For the investigation of other mesonic states we used a general two- 
fermion operator: 

0(t) = ^0xXx,iXx+ K i,t+Ti^(x,t;x + Ki,t + Ti), (34) 



X 



AX 



V^x — V l , VI — il ; ^1) T\ — 0, 1 . 

Charge conjugation eigenstates are constructed in the usual way by using a 
linear combination of (|34]) and its charge conjugate. The possible values for 
K\,T\ and are restricted by the requirement that the operator should also 
be an eigenstate of the symmetry transformations and J 1 . To reduce the 
complexity of the correlation functions and hence the computational effort 
we did not construct eigenstates of the inversion J? for all quantum number 
combinations. Nevertheless by considering the results of the transfer matrix 
formalism it is in all cases possible to determine the continuum quantum 
numbers of the measured energies (see section Q). 

Because of the complexity of the mass spectrum we calculated correlation 
matrices: 

Cij(t) = (0^)0,(0)) . (35) 

The operators Oi(t), i = 1 . . . r, differ only by different values of the smearing 
parameters. We use smeared sources with the smearing function 5 Xo (x) 



calculated by a Jacobi iteration [39]: 
"Sxo (x) := ^x,xo ) 

Sff(x) := r ^(sir 1) (x)+a{^- 1) (x + 2) + ^- 1) (x-2)}) . 

(36) 
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Parallel transporters have to be inserted into ( |36| ) so that the function 5 X0 (x) 
obtains the correct gauge transformation properties. After N iterations the 
normalized smearing function S' xo (x) with approximate Gaussian shape is 
defined by 



Sk (x) 



S- 



(N) 
xo 



is* 



^2 s L(x) s *o ( x ) 



1/2 



The calculation of this function only on even/odd sites is necessary for KS- 



fermions [40|. The use of smeared sources is most successful for the cal- 
culation of the pion mass. With N = 20 which is chosen throughout and 
a = 0.1 it is possible to suppress strongly other states in the CV-function. 
The radius r of the Gaussian curve S Xo (x) defined by 



J> 2 St (x)S X0 (x) 



(37) 



in that case is r ~ 2.5 . 

The different values of a for the operators in ( |35| ) are chosen below a = 
0.1. The one-particle energies are extracted from the generalized eigenvalues 
of the correlation matrix Cij(t). Because of the rich mass spectrum it was 
in most cases necessary to enlarge the correlation matrix. 



5.2 Finite size effects and dispersion relation 

Besides the knowledge of the mass spectrum a good control over the lattice 
artifacts is necessary as far as the elastic scattering phases are concerned. 
The finite size effects which are used to determine the scattering phases have 
to be separated from unwanted exponentially suppressed finite size effects 
(polarization effects). Furthermore, the correct dispersion relation has to 
be chosen to calculate the particle momenta of the scattering particles (see 



section 6A) and hence the scattering phases. 

For massive scalar bosonic quantum field theories Luscher derived a for- 
mula for the difference of a particle mass in finite (mi) and infinite (moo) 



volume [41, 42 1. For the lightest particle in the spectrum one obtains in two 



dimensions: 

Am = mi — m c 



-F(0) — — ==e- L *ff , L eff :=L moo . (38) 



TTlr. 



4 ^fhTL 



eff 
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The scattering amplitude F(y) depends on the variable v = (cj(p)Lj(q) — 
pq)/moo with p and q being the momenta of the incoming particles and 
u(p) the relativistic mass of the particles: 



w(p) = Vml, + P 2 • (39) 

Equation (^) is valid provided that no one-particle exchange scattering 
processes occur in the quantum field theory. This is fulfilled in the sine- 



Gordon model for (3gG = v27r. From the definition of F(v) in ref. |4l| one 
obtains with v = cosh(0) a connection between F[y) and the scattering 
matrix element S(6) in @. For v = one obtains F(0) tx —m^. Inserting 
this in @) one obtains 



Am = + A^^e- L *ff (40) 

with a numerical prefactor A. Considering the connection between the sine- 
Gordon model and the pion sector of the Schwinger model we use the formula 
([icf ) as a possible approximation to describe the mass shift of the pion in the 
Schwinger model due to finite size corrections. Because of additional states 
in the spectrum the formula ( f40"l) which we use for fits is not as theoretically 
justified as Liischer's formula (|38|). 

The numerical results for the finite size effects on the pion mass are 
shown in fig. || for j3 = 5 and (5 = 10. The value of the ratio mo/e = mo\[fi 
is fixed so that corrections to F(0) with respect to mo/e and therefore to the 
amplitude A are the same in both cases. A fit to the data in terms of the 
ansatz ( |40j ) shows a good agreement between this formula and the numerical 
data with similar amplitudes for both couplings (see fig. ||). It turns out 
that strong finite size effects exist for L e jf < 6. For £ e // > 8 the mass shifts 
are far below 1%. In fig. § they are comparable to the statistical errors for 
L e ff = 8. Therefore for all Monte-Carlo simulations we have chosen the 
spatial extension so that the effective length is L e ff > 8. 

For the determination of the particle momenta and the elastic scattering 
phases of the pions a correct dispersion relation has to be used to reduce 
0(a)-errors. The results for the energy-momentum relation E(p) of the 
pion in fig. [5] can be compared with the continuum and the lattice dispersion 
relations. As expected the bosonic dispersion relation 

2sinh-J = (^2sinh— J + p 2 , p:=2sin| (41) 

gives the best description for the numerical data. 
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Figure 4: The pion mass as a function of the spatial lattice extension for two 
different couplings f3 and fixed moy/P = 0.075\/l0- The mass in the infinite volume 
moo results from a fit with the function (|40| ) to the numerical data. 




Figure 5: Comparison of the different theoretical dispersion relations with the 
numerical results for the pion energy (left). Ratio of the pion energy to the bosonic 
lattice dispersion relation for two different lattices (right). 
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A better resolution of the difference between numerical results and the 
bosonic dispersion relation is obtained by taking the quotient of the simu- 
lation data and the bosonic lattice dispersion relation (second diagram in 
fig. H). Significant deviations from the theoretical curve which are above 
one percent occur for both lattices for momenta p > 0.5. This has to be 
compared with the simulation results concerning the scattering phases in 



section 6.5. The momenta of the scattering pions with masses m n < 0.4 
are in the region p = 0.1 .. . 0.2. Errors which occur by using the bosonic 
dispersion relation are therefore negligible. 

5.3 Light mesons 

In the strongly coupled Schwinger model the (approximate) scaling of the 
pion mass m n (mo) is predicted (see section |2.2| ). With regard to the corre- 
sponding formula (0) we check asymptotic scaling in the numerical simula- 
tions. Equation (|?|) is also used to determine the strong coupling region in 
which the predictions from the sine-Gordon model are reliable. 

First numerical results for the pion mass with reasonable agreement with 
© for m /e < 0.1 were presented in ref. j38|. The results from ref. [p5| , 



fig.l] for (3 = 4 show that a scaling of the pion mass according to 

^ocf^V, p = 0.689(10) (42) 
e V e / 

exists for mo/e up to 0.4. Nevertheless there are deviations for (3 = 2 
and (3 = 4 in ref. [25, fig.l] from the analytically predicted curve. As 



the continuum limit is reached for (3 — > oo these deviations decrease for 
higher values of the coupling (3 as expected. For = 10 and Q = the 
simulation data for the pion mass are compared with the analytical ones (||) 
and (^) in the double logarithmic plot in fig. |(| The data for mo/e < 0.25 
nearly coincide with the semiclassical formula (||) and are slightly above 
the exact prediction (|^). From these small deviations between numerical 
data and analytical formulae one can conclude that for couplings (3 > 10 
the expectation value on the lattice is dominated by the topological sector 
Q = 0. It is therefore possible to neglect non-trivial topological sectors for 
high values of f3. In contrast to the region mo/e < 0.25 it is visible in fig. ^ 
that for mo/e > 0.4 the pion mass does not scale very well according to ([?]). 

Besides the pion a scalar singlet with I PG = ++ ("/o" -meson) is pre- 
dicted in the sine-Gordon model for (3sg = V^/k. In the simulations a singlet 
with the correct quantum numbers can be identified. To compare the data 



with the predictions from section 2.2 we used the dimensionless mass ratio 
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Figure 6: Comparison of the numerical results for m^mo) {L x T = 48 x 48 to 
86 x 86) with the analytical predictions (||) and The parameters marked with 
dashed lines are the simulation points that have been used for the determination 



of the scattering phases (see section 6.5) 



mf /m 1T . This ratio is expected from eq. (||) to be mf^jm^ = \/3. The 
results for the /o and pion mass for [3 = 5 and (3 = 7 are plotted in fig. ^. 
For a bare mass of about mo — 0.14 the mass ratio from the simulations 
agrees with the predicted value. But for smaller values of the bare mass 
there are deviations from the \/3-line in fig. |?]. 

The third predicted particle in the Schwinger model is a pseudoscalar 
singlet ("77" -meson). In the massless case uiq = its mass is m^/e = \/2/tt 
(eq. (|8|)). It follows from dimensional arguments that corrections to this 
value in the massive model should depend on m-o/e only. Fixing the ratio 
mo/e in lattice simulations it turns out that for couplings (3 > 5 the ratio 
m v /e is indeed nearly constant (see fig. |8|). The simplest ansatz for the 
dependence of the 77-mass on the parameters in the Schwinger model is 
therefore: 



e V 7r v e / 



(43) 



This function can be used to fit the numerical results for m v for fixed cou- 
pling and various bare masses. The results for f3 = 5 and j3 = 7 in fig. show 
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Figure 7: Mass ratio mf /m w for (3 = 5.0 and [3 = 7.0. The expected value for the 
mass ratio for m /e — > according to eq. (^|) is mj /m T = \/3 (dashed lines). 
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Figure 8: Mass of the ?7-meson for two different operators in the same lattice sector. 
Three different values of (3 and mo are used with a fixed ratio mo/e ~ 0.16V7. 
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Figure 9: Dependence of the 77-mass on mo/e for (3 = 5 and /3 = 7. The ratio 
/ e(mo / e) according to formula (R) is plotted for comparison. 



that in both cases a nearly linear dependence with p ~ 1 and similar ampli- 
tudes A is obtained. Because of the value of p and A the 77-mass is larger 
than the mass of the pion in the parameter region considered (mo/e < 0.5). 
This statement is also correct for the simulation points used for the deter- 
mination of the scattering phases. Hence the lowest two-particle energies 
for these simulation points result from a two-pion system. 



Analogously to the pion (see section 3.1) we expect that the masses 
of the /o~ and 77-meson for different topological sectors are significantly 
different for large (3. Therefore the calculated masses for different ensembles 
with a different probability distribution p(Q) but same parameters might 
be different. We expect that the resulting errors for the data, especially 
for (3 = 1 and (3 = 9 in figs. [/],|| and ^, are small due to the fact that all 
simulations were started with Q = and hence the sector Q = strongly 
dominates in ensembles for (3 = 7 and = 9. 



5.4 Heavy Mesons 

In addition to the predicted particles in the Schwinger model we observe a 
large number of other scalar and pseudoscalar particles. The lowest states 



which could be determined with sufficient accuracy are shown in figs. 10 and 
|Tl| . The masses of these states are above the masses of the three predicted 
mesons. 

Many of these states might be unstable. For example, it follows from 
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Figure 10: Pseudoscalar triplets with a G = +1 (I PG = 1~+ = T~ + ) (left) and 
scalar triplets with a G = -1 (I PG = l +_ = T+~) (right). The dotted lines 
represent the sum of the masses of those particles which might be possible decay 
products of the plotted states. 



fig. |l(] that the Trf G = T 2 _+ -particle might decay into a w- and a Jo- 
meson. The states which are definitely stable in the examined parameter 
region are the Tf + - T+-- and the 5 ! -meson. That means that these 
particles do not have decay channels as far as the determined mesonic states 
are concerned. There might be other particles like triplets with quantum 
numbers (ap = +1, oq = +1) which have only a slight overlap with the 
operators we used in our simulations. 



In figs. 10 and 11 masses with m > 1 are to be understood mainly 



qualitatively, because one expects substantial 0(a)-corrections. 

Summarizing the results for the mass spectrum we see that we are able to 
identify the predicted mesons — ir, fo and r\ — and their masses in Monte- 
Carlo simulations. The analytical predictions for the scaling of the pion- and 
/o-mass (see section 22) are confirmed with small deviations. Best results 
are obtained for (5 = 10 whereas for (5 < 5 deviations from the continuum 
results occur. Considering the measurements for the various couplings the 
predictions from the sine-Gordon model for the strong coupling region of 
the Schwinger model are reliable for mo/e < 0.25. 

Concerning the 7T-7T scattering those particles with m ~ 2m n are of 
interest because they might occur as resonances in the elastic scattering 
processes or they might be confused with the two-pion-energies in finite 
volume as it will be shown in section 6.3. 
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Figure 11: Left: Different particles on the lattice which correspond to singlet states 
in the continuum (I PG = PG = S p ). Right: Complete overview of the measured 
states apart from the three light mesons n, f and r\. 



6 Scattering phases 

The aim of this section is to present the methods and results for the determi- 
nation of the scattering phases for the elastic 7r-7r scattering in the Schwinger 
model. As already mentioned in the introduction we use a formula which 
was derived by Liischer for massive quantum field theories |],||,|, 42]. This 



method is based on the fact that for not too small volumes the energies of the 
two-particle states in finite volume are determined by the elastic S-matrix 
of the model in infinite volume. Therefore Liischer's method is suitable for 
the investigation of scattering processes on the lattice with Monte-Carlo 
simulations. 



6.1 Liischer's method 

The system we are interested in consists of two bosonic particles of mass 
m in 1 + 1 dimensions in a finite box with spatial length L. We work in 
the center-of-mass system with zero total momentum^. The formula for 
the determination of the scattering phases 5 for the elastic scattering of the 
particles in infinite volume was given in ref. 0]: 

25(k) = -kL mod 2vr . (44) 

3 An extension of Liischer's method to systems with non-zero total momentum was 
proposed and successfully applied in ref. Jlo[. 
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k is the momentum of one particle in the finite box in the center-of-mass 
system. It is related to the discrete two-particle energy E by the relativistic 
formula: 



E = 2\/m 2 + k 2 . (45) 

For a potential with short range it can be easily shown that equation ( |4"4| ) 
is exact in quantum mechanics Q. Also for a bosonic scalar quantum field 
theory in four dimensions a relation between the scattering phases and the 
particle momenta in finite volume can be derived Q. 

An exact derivation of (^) in an arbitrary quantum field theory is in 
general not possible. Nevertheless it is plausible that fl44|) holds provided 
that § 

• the examined quantum field theory has no massless particles in its 
spectrum 

• the mixing of the relevant two-particle states with other many-particle 
states is excluded 



the range of the two-particle interaction is much smaller than half the 
spatial extension. 



In the sine-Gordon model no inelastic scattering processes exist [j!5[|. Never- 
theless such inelastic processes are possible in the Schwinger model. There- 
fore in the Schwinger model one has to ensure that the lowest two-particle 
energies are below the inelastic threshold E = Am. Furthermore, one has to 
guarantee that no production of other particles occurs in the relevant energy 
region. 

If the extension of the spatial volume is too small, further exponentially 
decaying finite size effects (polarization effects) appear in the two-particle 



energies. This leads to corrections to formula (44). As far as the scattering 
phases are concerned it is difficult to determine numerically these finite size 
effects. To reduce these unwanted effects on the two-particle energies the 
minimum lattice size for all simulations is chosen so that the finite size effect 
for the pion mass is negligible. 



6.2 Correlation functions 

For the determination of the two-particle energies in the finite volume we 
start with the successful ansatz of refs. P,|9|)P!2^|. This means that we con- 
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struct a correlation matrix with time-slice-operators where different opera- 
tors are characterized by different orthogonal functions: 



Of TT {t) = ^2</> x ,yXx,tXx,t Xy,t+TXy,t+TUi(y ~ x) , (46) 



L 2 



6 x , y = (±l) x (±l) y , r = 0,l, 
cos(pjx) 



Wj(x) 



sin(pjx) , pi = ^-i, i = 0,l, 



These operators having the property that the fermion and antifermion in 
each meson operator are local in time are denoted by LT = LT$ (r = 0) and 
LT\ (r = +1). The definition of the functions Wj(x) and 4>^ :Y depends on 
the choice of the lattice symmetry sector. In some sectors it is also possible 
to use two definitions for ^ x>y : </>x,y = (— °"i) x ( — l) y or <^ x>y = cri x . The 
corresponding operators are denoted with a 1 and 2 respectively appended 
to the quantum numbers. 

The operator ( f46| ) has the quantum number oc = +1- To construct 
eigenstates on the lattice which are odd under charge conjugation we use 

Of LT {t) = J2J2 ^ Xx,t+lXx,* Xy,tXy,t+i ^(x, t)U* 2 {y, t)w f (y - x) , (47) 
x,y 

which in general is not an eigenstate under charge conjugation. This is an 
operator where the fermion and antifermion in each meson are non-local in 
time (NLT). Like in the case of the LT-operators the definition of Wj(y — x) 
and </> Xj y is determined by the symmetry sector. 

The advantage of the special construction (^7|) is visible if one constructs 
an operator with a definite quantum number ac' 

= Of LT (t) + a c (%Of LT )(t) 

= ^Xx,t+iXx,tXy,tXy,t+i U 2 {x,t)U2{y,t)uJi(y - x) 
x,y 

• {</>x,y + 0-C07</>y,x} . (48) 



Because of the expression (|48| ) the calculation of Green-functions with the 
TVXT-operator is not more expensive than the calculation with the LT— 
operator. 

As it will be shown in the following section the four-meson correlation 
functions of the operators ( |46"| ) and ([47]) include contributions from one- 
particle states. To avoid such contributions for the determination of the 
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7r-7r scattering phases it appears advisable to calculate correlation functions 
which are constructed according to the (/ = 2)-sector in the continuum 
because there no one-particle states occur. Details will be seen in section 



5.3| . The procedure of the derivation of these correlation functions is the 
same as for the CV-function: The starting-point here are the expectation 
values C con t,ij(t) = (Oi(t)Oj(t)) in the continuum. The operators Oj and 
Oj are defined by 

Oi(t) = ■^y'dxdy^(y-x)V; 1 (x,t)75V' 2 (x,t)^ 1 (y,i)75V' 2 (y,t), 
Oj(t) = 7 ^y'(ixdy^-(y-x)V; 2 (x,t)75^ 1 (x,t)^ 2 (y,t)75V' 1 (y,t) • 

By naive translation of the correlation functions C con t t ij(t) into the lattice 
formulation one obtains among others the following terms 

x,y,z 

-2Re [M*- 1 (1,4)M- 1 (1,S)M*- 1 {2,3)M- 1 (2,4)] ) u 
•(-l) 1+2+3+ V(y - x)^-(z) x , y z , o , (49) 

l = (x,t), 2 = (y,t + r), 
3 = (z,0), 4 = (0,r'). 



The prefactors in (|49| ) are chosen in such a way that the expression (^) coin- 
cides with corresponding terms in the correlation function of the operators 
(p6|). In the following the index W7r always denotes a correlation function 
of the type (|49|). For r = r' = 0, 1 the expression ( |4^ ) is related to the 
LT, LTi-operator and for r = 1, r' = to an LT\-LT correlation function. 
If one takes the naive continuum limit of these lattice correlation functions 
it turns out that only the expression ( (49| ) contributes to the expectation val- 
ues C con t t ij(t). Hence we expect the C^n (t)-function to yield the dominant 
contribution to the 1 = 2 continuum sector and therefore to be sufficient to 
calculate the two-particle energies in this sector. 

The advantage of the expression (H) is obvious. Compared to the calcu- 
lation of e.g. (O lt (t)0 LT (0)) the calculation of ( |49j ) needs merely a fraction 
of the conjugate gradient calculations. 
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Figure 12: Results for the pion mass for different four-meson correlation functions 
in different lattice sectors. The dashed line is a fit to the data, whereas the solid 
line represents the mass obtained from the two-pion correlation function C n . 



6.3 One particle states in four meson correlation functions 

The numerical results for the four-meson correlation functions with the oper- 
ators ( |46| ) and (|47]) show that in nearly all sectors the main contributions and 
hence the lowest energies can be identified with the single-particle masses of 
the model. Just as for the two-meson correlation functions the pion-state 
here also occurs with the best signal-to-noise ratio. The comparison of the 
determined masses with the pion-mass calculated with the two-pion corre- 
lation function C v (t) shows a good correspondence between the results for 
the two-meson and four-meson correlation functions (see fig. [12] ). 

This result can be explained in terms of the general properties of the 
SU(2) flavour group of the Schwinger model: The structure of the group 
SU(2) is such that some irreducible representations may be expressed either 
in terms of tensors of rank four or tensors of rank two. Therefore it is 
possible to examine properties of particles with isospin 1, i.e. triplets, with 
fermionic eight-point functions. According to table [l] the lattice sectors 
which are connected to a I PG = 1 h continuum state are (a\,aj,ac) = 

( h), ( — h+) and (+ + +). In fact we observe the pion in all these sectors 

(see fig. |l^). Besides the pion we have determined also other mesons — the 
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fo— and T x -meson — with appropriate fermionic eight -point functions. 

According to table |l] the pion as well as the tt-tt states occur in most 
sectors with ac = +1. Therefore the lowest energy and consequently (pre- 
sumably) also the signal-to-noise ratio in these sectors is dominated by the 
pion. Furthermore, according to the results of section [TJ there exist pseu- 
doscalar mesons with masses m ~ 2m n . These particles possibly could make 
the determination of the two-particle energies more difficult. To avoid these 
problems we used for the reasons given above the C 7r7r (t)-function ( |49| ) for 
the determination of the scattering phases. In fact the numerical results 
discussed in section |6.5| show that by using the CV^^-function the contri- 
butions of the single-particle states are suppressed. In particular the lowest 
energies in the C nn (i)-function are the tt-tt energies. 



6.4 Connected correlation functions 

One might suspect that the only cause for the appearance of the one-particle 
energies in some four-meson correlation functions are poles in the discon- 
nected parts of the correlation functions. To refute this hypothesis we exam- 
ined a four-meson correlation function C y (t) = / 0\ T (t)Of T (0) V (The fol- 
lowing discussion also holds for other operators (LT±, NLT, . . .) with more 
general space-time arguments of the fermions.) This correlation function 
separates into a sum of connected correlation functions and disconnected 
parts: 

Cij(t) = + disconnected graphs . (50) 

The shaded symbol represents the connected Green-function ((. . .) c ). The 
external legs denote quark-antiquark pairs in which the incoming mesons 
on the left live on the time slice t and the outgoing mesons on the right on 
the time slice 0. 

The disconnected parts consist of 




The dots in (^2[) represent graphs with three external legs like 

In the usual procedure only the vacuum parts (A) are subtracted from 
the Green-function to obtain the connected Green-function 

(O i (t)O i (0)) c = ; <O i (t)O i (0)>-(O i (t))<O J -(0)> . (53) 

In our simulations it turns out that the expressions (B) and (C) are not 
negligible and must also be subtracted in eq. (|53|). However it can be 
derived easily that (B) as well as (C) disappear for i ^ j because of the 
orthogonality of the functions Wj(y — x). For example, the term (B): 

(B) oc ^2 ( / , x,yW i (y-x)tJ i (y / -x / ) 

x,y,x',y' 

• < Xx,tXx,i Xx'A,Xx',t } < Xy,tXy,t Xy',t Xy',t } (54) 
can be rewritten to obtain 

L 



(B) oc 8ij— 0i,x cos(pix) ( Xx,tXx,i Xo,t Xo,t ) 

x 

• ^2 COS (P«y) ( Xy,tXy,t X0,t„X0,t o ) > (55) 

y 

</>x,y = :01,x02,y > * 0, ~Z ■ 

In the same way one obtains (C) oc 5ij. This means that the disconnected 
parts (B) and (C) contribute to the diagonal elements of the correlation 
matrix only. 

This result shows that the one-particle contributions to the four-meson 
correlation functions can not come from the terms in (pT]) only, because we 
obtain in our simulations that the pion state also occurs in the non-diagonal 
elements of the correlation matrices (see table ^) . 

Actually, the disconnected parts are a nuisance for the determination of 



two-meson energies: From equation (55) it follows that the disconnected 



parts are products of propagators of moving mesons. Therefore from 
one obtains energies which are twice the energy of a single free meson with 
momentum pj. In fact we could determine these energies with a good signal 
in the Green-function Cij(t). The overall results for Cij(t) as well as the 
measured energies from the disconnected parts (B) + (C) alone are listed in 
table ^. Comparing the energies which correspond to each other from both 
measurements good agreement within the errors is visible. 
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/ 0.503(14) 0.506(6) 0.492(6) 0.498(17) \ 

0.505(4) 0.50(2) 0.504(12) 0.516(6) 

0.498(2) 0.50(4) 0.495(9) 0.523(9) 

\ 0.494(2) 0.51(4) 0.495(9) 0.526(9) / 

Table 2: Pion-masses which result from each of the 16 matrix elements for a 4 x 4 
correlation matrix. The values are obtained with an LT operator for [3 = 5.0, 

m = 0.16 and L x T = 24 x 64 in the sector {<j\, 07, oc) = ( For 

comparison: The pion-mass for a single-particle operator is m n = 0.4991(3). 

These energies from the disconnected parts make it more difficult to 
analyse the simulation data and to determine two-meson energies we are 
mainly interested in. In some cases they might even be confused with the 
two-pion energies. In order to avoid contributions from the expressions (B) 
and (C) we considered only correlations with i ^ j, because as we have 
seen above all non-diagonal elements of (B) and (C) vanish. In fact we 
calculated r x r correlation matrices D(t) which are not symmetric in their 
indices: 

Dij(t) := Cj iT .+j(t) 

= {Oi(t)O r+j (0)) , i,j = l...r. (56) 

Moreover, we did not use momenta pj = 0, it for the wave functions. In this 
case it can be shown that the contributions from (^) vanish, too. 

Altogether we use the function D^it) for the calculation of the two- 
pion energies. This is a correlation matrix (|56|), for which according to (|49|), 
only those parts contribute which are relevant for the 11-11 scattering in the 
(/ = 2)-sector. 

Analyzing the numerical data for D n7T (t) it turns out that the best re- 
sults are not obtained by calculating the generalized eigenvalues but just by 
diagonalizing the correlation matrix. In this case the eigenvalues show the 
following behaviour (T — > 00): 

A,(t) * =°° diatf e~ tE > [1 + 0(e- tAE <)] , a t = ±1 . (57) 

Like in eq. (|l9|) AEi denotes the smallest difference between E\ and other 
spectral energies. Following closely the proof in ref. Q, equation fl57| ) can 
be derived easily for the D n7T (t) -cor relation matrix. 
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Cij{t) 


(B) + (C) 


2£ , feos (m 7r ,p i ) 




0.506(5) 






Energy E2 


1.10(3) 


1.21(4) 


1.124(9) 


Energy E 3 


1.53(11) 


1.49(3) 


1.418(7) 


Energy E± 


1.89(16) 


1.89(4) 


1.776(5) 



Table 3: Comparison of the results from the four-meson correlation matrix Cy(t) 
with the values resulting from the disconnected parts (B) + (C). The energies were 
obtained for the same simulation parameters (3 = 5.0, mo = 0.16, L x T — 24 x 64 
from a 4x4 correlation matrix. Additionally the energies which were calculated with 
the bosonic dispersion relation are listed (m, — 0.506(5), = 2iri/ L 1 i = 1 ... 3). 



6.5 Results for the scattering phases 

For the determination of the two-particle energies we calculated a 4 x 4 
D niT (t)-matrix with the LT-, LT\ -operators (see section |6.2| ) and the LT\— 
LT operator combination for all relevant lattice sectors. The simulation 
results show that a clear signal for the 7r-7r energies occurs only for the LT- 
operator in the (01, 07, ac) = (+ + +)2 lattice sector. For the calculation 
of the energies in this case it was frequently necessary to enlarge the corre- 
lation matrix to obtain a better separation of the different energies. Hence 
we extracted the lowest energy only whereas higher energies could not be 
determined with sufficient accuracy. 

To improve the signal for the lowest energy we use smeared operators 
(see section |5.l| ). The optimal value for the smearing parameter a from 
simulations with low statistics turned out to be a ~ 0.02. Here one has 
to be aware of the fact that the replacement of the antifermions in a two- 
meson operator 0{t) with quantum number ac(0) by smeared antifermions 
("smeared source") leads to a smeared operator O (t) which has no definite 
behaviour under charge conjugation. In such cases it is useful to chose a so 
small to have dominant contributions from the sector with ac = ac{0). 

The results for the lowest eigenvalues of two typical correlation matrices 
are shown in fig. [l^. Sufficiently good plateaus in the effective energies are 
obtained for relatively high values of the time argument t. For small values 
of t the eigenvalues show a qualitatively different behaviour. Therefore the 
calculation of generalized eigenvalues is not successful here, because for the 
method of generalized eigenvalues the reference time to of equation ( Jl8| ) has 
to be chosen small to reduce the statistical errors. 
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Figure 13: Two-particle energies for ((3 — 10, mo = 0.11 ) for two different spatial 
extensions of the lattice in the sector (<7i, <Ti, ac) = (+ ++)2- The lowest eigenvalue 
of the correlation matrix and its effective energy are plotted. 



For the -D^^-function the use of smeared sources is not as advanta- 
geous as for the two-pion correlation function C n (t). Despite smearing the 
matrix elements of the subtracted function D^ 7r (t):=D^'J(t) — D^J(t + 2) 
for t = 10 and (/? = 10, too = 0.11) have a bad signal-to-noise ratio of 
about 0.02 to 0.08. The subtracted correlation function is used here in or- 
der to avoid any contributions from constant terms. The signal-to-noise 
ratio above has to be compared with the signal-to-noise ratio of 4.18(2) for 
the subtracted function Cs^i^C^t) — C w (t + 2) with a = 0.1. Therefore 
we had to generate for each point a minimum of 65000 independent config- 
urations to reduce the statistical errors of the two-pion energies to about 
1%. 

The comparison of the -D^(i)-function with the correlation function 
C LT (t) = (Of T (t)Of T (0)) justifies the usage of the former function: The 
ratio D^ i j(t)/Cj< T (t) of the matrix elements of both functions for the pa- 
rameters given above is less than 10~ 3 . This means that the D^(t) parts 
of the C LT (t)-correlation function are strongly suppressed. Thus it is ap- 
propriate to use the correlations D^(t) from the beginning. 

The scattering phases which were calculated from the two-particle en- 
ergies by the bosonic dispersion relation ( |4l"| ) are shown in figs. [l4| and |l5]. 
In fig. [14| the numerical values confirm well within errors the analytical pre- 
diction (g) for the 7T-7T scattering. This means that in both cases (/3 = 10, 
too = 0.075) and ((3 = 10, too = 0.11) for a ratio too / c = moy/P of 0.24 
and 0.35 respectively the strong-coupling prediction describes the data very 
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Figure 14: Scattering phases for two different bare masses and different lattice 
sizes. The analytical prediction for the 7r-7r scattering according to eq. (^) is 
represented by the solid curve. The mass of the scattered pions is = 0.2545(10) 
and m-n — 0.3367(6) respectively. The solid/dotted bars in the lower diagram are 
the statistical errors of the data for T — 48 /T — 64. 
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Figure 15: Scattering phases for j3 = 10 and too = 0.14. The pion mass is 
= 0.4010(2). 

well. Therefore the analytical strong-coupling prediction for the elastic tt-tt 
scattering phases in the Schwinger model as resulting from the sine-Gordon 
model are confirmed. 

If the ratio mo/e is increased further to mo/e = 0.44 for {(3 = 10, tuq = 
0.14) it becomes obvious from fig. ^ that the prediction (g) is no longer 
valid. Indeed the lowest numerical determined scattering phases in fig. ^ 
lie below the analytical curve. Considering our numerical results for the 
mass spectrum these deviations can not result from inelastic processes. The 
threshold for possible inelastic processes like tttt — > 7777 or nir — > f fo for the 
parameters in question is far above the determined two-particle energies. 
Furthermore, we did not find any particles, e.g. (I PG = ++ )-mesons, 
which might occur as resonances in the examined energy region and the 
lattice sectors considered. Therefore the deviations in fig. |l5| can not be 
caused by resonances in the elastic scattering channel. 

Because there are no inelastic thresholds near or below respectively the 
calculated energies we can conclude that the investigated pion-pion state is 
a definite two-particle state which does not mix with other many-particle 
states. Therefore according to section ^T] the applicability of Liischer's 
method is guaranteed. Within the statistical errors the numerical results do 
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not indicate that it is necessary to consider any CDD-factors for the elastic 
pion-pion scattering. Nevertheless it is possible to multiply the "minimal" 
solution (JoJ) with CDD-factors of the type (|TT[) , e.g. two factors with a\ — 
—ct2, without changing the shape of the scattering phases qualitatively. 



7 Summary 

In this paper we present the mass spectrum of light mesons in the massive 
Schwinger model with an SU(2)j flavour symmetry. Furthermore, Liischer's 
method is successfully applied to the calculation of mesonic scattering phases 
in this model. For strong coupling the analytical predictions for the pion 
mass and for the elastic scattering phases of the pion-pion scattering are 
confirmed. 

Investigating the mass spectrum in the Schwinger model it is shown that 
for (3 = 10 and for small bare fermion masses mo the dependence of the pion 
mass on the bare parameters corresponds to the analytical predictions. A 
significant improvement of the data is achieved compared to the results for 
/3 = 4 [25|. This was expected as the continuum limit requires j3 — > oo. 



The /o and r\ singlet states which are analytically predicted in addition 
to the pion are investigated in the strong coupling region, i.e. for small 
values of m,Q^[]5. The numerical results for the mass of the /o-meson are 
for (5 = 5 and (3 = 7 in the region of the analytically predicted value of 
rrif = V3m n . The mass of the 77-meson is predicted to be = 
in the massless SU(2)j Schwinger model (e is the coupling parameter in the 
continuum). We find in the massive case that there are corrections to this 
value which are linear in m,Q. 

Additionally in the mass region [m n ,4:m n ] we find a rich mass spectrum 
consisting of scalar and pseudoscalar triplet and singlet states. We assume 
the lightest of these particles to be stable. 

The results for the elastic scattering phases of the pion-pion scattering 
show that the relevant two-pion states have only a slight overlap with the 
operators considered. In fact the signal-to-noise ratio of the correlation 
function used is below 10% despite smearing. To improve the separation of 
the pion-pion energies from other states with large amplitude we use im- 
proved techniques like specific correlation functions for the two-pion states 
and enlarged correlation matrices. In this way it is possible to obtain scatter- 
ing phases for the elastic pion-pion scattering using the method of Liischer. 
For strong coupling these scattering phases are in good agreement with the 
analytical predictions. 

The experience gained in this project may be useful for the calculation of 
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elastic scattering phases in other gauge field theories. It turns out that the 
investigation of scattering processes with composite particles like mesonic 
states is much more difficult than the scattering of particles which are no 
bound states. This is in agreement with investigations of other fermionic 
models pLpLpj. The main problem is the poor overlap of the two-meson 
states with the conventional operators. To compensate for this we had 
to generate up to 100000 independent configurations per simulation point. 
Therefore the development of improved operators should be the main task for 
investigations of mesonic scattering processes in QCD by Luscher's method, 
e.g. of the /^-resonance in a tt-tt scattering. 
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